Universal spatial correlations in random spinor fields 
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We identify universal spatial fluctuations in systems with non trivial spin dynamics. To this end 
we calculate by exact numerical diagonalization a variety of experimentally relevant correlations 
between spinor amplitudes, spin polarizations and spin currents both in the bulk and near the 
boundary of a confined two-dimensional clean electron gas in the presence of spin-orbit interaction 
and a single magnetic impurity. We support or claim of universality with the excellent agreement 
between the numerical results and system-independent spatial correlations of a random field defined 
on both the spatial and spin degrees of freedom. A rigorous identity relating our universal predictions 
with response functions provides a direct physical interpretation of our results in the framework of 
linear response theory. 
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The field of spintronics, which deals with the use of the 
electron spin of freedom as means to transmit, store and 
process energy and information, has experienced impres- 
sive advances during the last decade The possibility 
of manipulating spin densities and currents by means of 
their correlations with the easier to control charge (and 
electrical current) degrees of freedom has been in the fo- 
cus of semiconductor-based spintronics. Here the spin 
Hall effect [2], the creation of a spin imbalance across a 
sample generated by a charge current in the presence of 
spin-orbit interaction (SOI), is a primary example of such 
spatial correlation and the close connection between the 
spin and charge degrees of freedom. 

In this context, universality of spatial correlations can 
arise from two basic mechanisms. In systems where the 
mean free path If is much smaller than the system size 
L, average over the random distribution of obstacles pro- 
duces results which are largely independent of the geom- 
etry of the confinement. This diffusive limit has been 
extensively studied using diagramatic techniques based 
on disorder average |4[ which, however, cannot deal with 
ballistic systems where If /L ^> 1, a regime which is now 
easily achieved in high-mobility semiconductor 2D elec- 
tron gases [3]. In this ballistic case universality arises 
due to electron scattering with the irregular boundaries, 
namely, from the presence of classical chaos. It is in this 
regime where semiclassical approaches to universality [5] 
(and its breakdown |6[) in spin and charge transport in 
the presence of SOI have been very succesful. 

Spatial fluctuations in ballistic spinor systems have 
been addressed in [7( in the limit of vanishing SOI, and 
recently in [8[ the spatial correlations of charge densi- 
ties for the bulk have been studied by means of Random 
Matrix Theory. The extension to systems with non-zero 
local spin polarization (which is realized by an STM tip, 
for example) and in the presence of boundaries requires 
substantial technical and conceptual steps beyond Refs. 



[7, 8]. Our goal is to fill this gap. 

Our starting point is the Schrodinger equation (1(1) is 
the unit operator in spin (position) space) 
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with boundary conditions 
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where we use ^ n (r, s) = ((r| ® (s|)|^ n ) and 5 S , S > is the 
kronecker delta. Eqs. (jl|2p describe the dynamics of an 
electron with effective mass m inside a quantum dot with 
Dirichlet boundary conditions along in the presence of 
SOI with inverse precession length fc Si0 , and of a magnetic 
impurity (polarized in the \s =t) direction) located at po- 
sition q. The momentum operator p = —ih(d/dx,d/dy) 
acts on the orbital degrees of freedom r = (x,y), and 
(d- x ,a y ,& z ) are the Pauli matrices acting on the space 
spaned by the eigenstates of a z (denoted by \s) with 

*=t,l). 

Time Reversal Invariance (TRI) is expressed by the 
condition [9] HT = TH, where T = —icr y K is the time 
reversal operator and K indicates complex conjugation 
in the eigenbasis of o z . Since Dirichlet boundary con- 
ditions are also invariant under time reversal, without a 
magnetic impurity eigenstates of H come in degenerate 
(Kramers) pairs (|^ n ), |^n) = T\ty n )). Kramers degen- 
eracy prevents the very existence of non-zero local spin 
polarization if the state of the system is given by an in- 
coherent (statistical) superposition, 

Tr pn&j = for p n ~ |*„)<*„| + (3) 

This is the reason why Ref. [8[ deals only with correla- 
tions of the local charge, as this is the only density that 
does not vanish when averaged over the Kramers pair. 
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However, any specific coherent superposition a\^ n ) + 
b\ty^) shows indeed a very rich spatial and spin structure, 
which is washed out if we neglect coherences. Physical 
realizations of such coherent superpositions are achieved 
by breaking TRI, and in particular by introducing spin 
selective boundary conditions, as in Eq. (j2j). As an exam- 
ple, Fig. (pQ) shows the spin polarization in z-direction for 
a typical eigenstate of our hamiltonian where the mag- 
netic impurity is localized inside of a ballistic cavity. Here 
we will present a theoretical approach to understand and 
predict the spatial statistics of such imprinted pattern. 
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FIG. 1: Spin polarization in z-direction in units of L~ 2 for 
the linear combination of the two degenerate states n = 400 
of the " star billiard" such that the combined state is polarized 
in the | t) direction at the origin (k s . Q L = 10). 

For the numerical simulations we consider the desym- 
metrized "star billiard" (shown in Fig. [I]), which is char- 
acterized by only one parameter (the radius of the cir- 
cular arc in units of the length L of the horizontal seg- 
ment). For /c s .o = this billiard is known to display 
hard chaos in the classical limit and therefore it is ex- 
pected to exhibit universality in the spatial correlations 
of its quantum eigenstates [11[. In order to construct 
local observables, we explicitely diagonalize the hamil- 
tonian using the method reported in [12], which gives 
the corresponding Kramers pair for each eigenenergy E n 
. For a given position q of the spin impurity the (now 
unique) eigenstate of the system is given by the linear 
combination 



|C' b )=«(q)l*n)+6(q)|C), 



(4) 



where the coefficients a(q), 6(q) must fulfil |a| 2 + \b\ 2 = 1 
and are adjusted such that \£ n (q, s) = S s ^. For fixed n 
any local function ,F(^J' b (r, s), ^J' b (r / , s')) of the state 
will fluctuate when q is randomly choosen inside the bil- 
liard. This quasi-random character of the spatial fluctu- 
ations is used to replace the average over the impurity 
position q by an average over (a, b) on the unit sphere, 



^(C (q) ' b(q) (r,s),C (q),Kq) (r , ,s , ))dq 

= / 

J\a\ 



(5) 



'|a| 2 + |6| 2 = l 

This is the way we use the numerical Kramers pair 
(l^n)? l^n)) t° construct the expectation values of lo- 



cal observables T at fixed energy. Now we proceed in 
the usual way one studies spatial fluctuations of wave- 
functions in classically chaotic systems. We use the ex- 
act numerical results for local observables to perform an 
energy average where we expect universality to emerge. 
That this is actually the case can be seen in Fig. [2j where 
the two-point spatial correlation (^' 6 (r, s)ipn ,b (r\ s')*) 
of the spinor amplitude (calculated by impurity and en- 
ergy average of the exact eigenfunctions) is compared 
with the universal prediction of the Spinor Random Wave 
Model (SRWM) to be presented bellow. 
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FIG. 2: Comparison between the exact (left panel) spatial 
two-point correlation function (^(r, t)'0n(r / , t)*) (denned by 
impurity and energy average over 300 states near n — 3000) 
as a function of r, near a wall at x = (with r' = (0, 0.05L)) , 



see text. We use k s . /k 



: 0.1 and kL ■ 



and the universal results Cq 1 (r, r'; E) + Cj' (r, r ; ; E) (right), 



71* ' 
300. 



Any theory which attempts to describe the universal- 
ity of spatial correlations must deal with two different 
aspects: first, the issue of the universal behaviour of the 
amplitude correlator both in the bulk and near a bound- 
ary, and second the appropiate description of correlators 
beyond the bilineal form in the amplitudes. The origi- 
nal approaches to universal spatial fluctuations in chaotic 
systems considered these two phenomena to share the 
same origin (as both can be derived from Berry's ansatz 
stating that wavefunctions of classically chaotic systems 
behave as random superpositions of plane waves fill]). 
Later it was recognized [13] that the two-point correlator 
is actually an intrinsically microscopic object that can be 
derived without any further assumption from the exact 
Green function ^(E) = (H — £?±i0 + ) -1 by means of 
the formula 



1 1 

2ni q(E) 



<G-(r',r;£)-G+(r, r';£)>. (6) 



Here (...) denotes impurity and energy average (over a 
small window around E) and g(E) is the smooth part of 
the density of states. For systems with SOI and away 
from the boundary, the exact Green function is approx- 
imated by its bulk value to get the universal prediction 
for the bulk, 



C D (d;E) 
e ie C s (d;E) 



e-* 9 C s (d;E) 
C D (d;E) 



(7) 
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where C D (d;E) = fc + J (fc+d) + k_J (k_d), C s (d;E) = 
k+Ji(k+d) — k-Ji(k-d) and J n (x) are Bessel functions. 
We further defined k± = ^k 2 + fc| >0 ± fc s>0 with fc = 
^2mE/h 2 and r - r' = d(cos 0, sin (9) . The result can 
be obtained using a modified Berry ansatz including SOI, 
as in Refs. [!, LL41 ]. 

To go beyond the results for the bulk we use a mul- 
tiple reflection expansion to construct the matrix- valued 
Green function near a wall, assumed to be an infinite 
straight line at x = 0. Translational invariance in the 
x-direction suggests to perform a Fourier transform (in- 
dicated by a tilde) from x — x' to k x . The method pro- 
vides a closed form for the boundary contribution as (for 
notational convenience the dependence with E is kept 
implicit) 

Gi(k x ,y,y') = 2dG (k x ,y)(l - 2dG (k x ,0))- 1 G (k x ,y'), 

in terms of the bulk Green function 

, v-^ e -a±\y-y'\ , A x 
G (k x ,y-y') = 2^(±) — ^ [f±I - 9° y + ha x ) , 



and its normal derivatives dGo at the boundary. We 
introduced f± = g(k 2 — k±)/k X: and h = iga±sgn(y — y') 



with g = k x /(2^/k 2 + k 2 Q ) and a± = yj k 2 — k± (5fta± > 
0). 

In order to construct the correlation function in real 
space, the inverse Fourier transform of G\ is calculated 
in stationary phase approximation, well justified in the 
regime kL ^> 1. We obtain a linear combination of ex- 
pressions with phases $(k x ) = — k 2 • y + yjk \ — k 2 ■ 

y' + k x (x — x'), showing that the geometry of the sad- 
dle points $>'(k x ) = is a deformed version of the 
Snell law, with two (instead of just one) possibilities 
for the incoming and outgoing wavevectors. The limit 
k s .o/k —> 0, k s . Q L — >> const, where the SOI is not consid- 
ered for the stationary phase condition gives insufficient 
results for the spatial correlations, and we use instead 
a consistent solution for each independent combination 
of beams up to second order in k SmQ /k. The effect on 
the spatial correlations is then incorporated as a contri- 
bution Ci on top of the bulk result Co, Eq. (0. The 
effect of a nearby boundary on the two-point amplitude 
correlator is depicted in Fig. [2] showing excellent agree- 
ment between numerical simulations and our analytical 
(but very lengthy) formula. This boundary effect is the 
spatial analogue of the perimeter term in the Weyl for- 
mula [15] for the smooth part of the density of states of 
a system with SOI [l6| . 

We have checked that changes on the size of the energy 
window, relative positions r, r', position of the wall, SOI 
strength do not affect the quality of our results for all 
the entries of the correlation matrix. Therefore we con- 
clude that spatial correlations of spinor amplitudes are 



described by the formula ((6]), and universality emerges 
when the Green function can be approximated by its uni- 
versal limit for the bulk or near a hard wall. That this is 
precisely the case for chaotic quantum systems was shown 
for the scalar case in |17j . and the same argument (that 
paths with multiple reflections produce sub-dominant ef- 
fects) holds here. 

We now turn our atention to the experimentally more 
relevant case of spatial correlations for local densities. 
We consider observables of the form 



A d (r) = 5(r - r) <g> A. 



(8) 



For given position r, the choice A = 1 describes the local 
charge density, while A = &i gives the spin density in 
i-th direction. Introducing the spinor ip® ,b (r) = (r|^' 6 ), 
the numerical spatial density correlations are then con- 
structed from the numerical eigenstates by impurity and 
energy average 



C AB (r, r';£) = 



(9) 



Contrary to the amplitude correlators, density correla- 
tions are not bilineal in the components of the state and 
therefore cannot be directly related with the Green func- 
tion. Following a well established procedure in systems 
without spin, we assume that the spinor amplitudes have 
Gaussian fluctuations [13], and we replace the energy and 
impurity averages by a single average over a functional 
distribution of spinor fields, 



(r,r / ;£;)= / P(^)^(r) t i^(r)^(r / ) t ^(r / )^M. 



n SRWM ( 



The probability distribution P(ip) is Gaussian and there- 
fore uniquely given by its two-point correlation function 
/ P(^)^(r)^(r)tX>[^] = C(r,r';£), which we replace by 
the universal amplitude correlator based on the micro- 
scopic Green function. 

Having at hand a Gaussian theory with known two- 
point correlators, we can decouple averages over higher 
order functionals of the state by straightforward use of 
Wick's theorem. For the particular case of local observ- 
ables, this easily gives [16| 

C^f M (r, r'; E) = Tr [iC(r, r'; E)BC{v\ r; E) 

+Tr [iC(r, r 7 ; E)] Tr \bC(t\ r; E)] (10) 

where the trace is over spin variables only. Eq. (flQ]) is 
the main result of this paper. It allows us to derive uni- 
versal results for spatial correlations of local observables 
after inserting the expressions Co,Ci of the correlator 
C(r, r'; E) for the bulk or its modification near the wall. 
In order to check the underlying Gaussian assumption, in 
Figs. I3|4I we compare the result ([TO]) with numerical re- 
sults based on Eq. (j9]). We find considerable agreement 
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even for the subtle patterns emerging from interference 
effects due to the boundary. 

The physical relevance of Cab can be extended be- 
yond its statistical interpretation by noticing that the 
connected part of C^ WM at the Fermi energy E = Ep 
can be rigorously related through 



lim 



<t>AB(t)- 



iwt 



-dt 



w 



(11) 



iTi AC(r,r';E F )BC(r',r;E F ) 



to the d.c. component of the dynamical response function 
<t>AB{t — t') ([A d (r,t), B d (r,t')]) describing the change 
of the expectation value of A d at time t when an infinites- 
imal perturbation affects the system at time t' through a 
coupling with the observable B d [To[ . Eq. ([TT]) relates an 
experimentally accessible quantity, the response function, 
with the correlator quantifying the statistical fluctuations 
of the random spinor field. 




FIG. 3: Spatial correlation C ay(Tz (r, r'; E) of spin densities 
in the bulk (top) and near a wall at x = (bottom) as 
function of r. Left panels: results of Eq. (|9]) based on nu- 
merically obtained eigenstates. Right panels: universal pre- 
diction C%™ m {y,y'-E), Eq. (TO]). We use r' = (0,0.05L), 
ks.o/k = 0.1 and kL = 300. Averages are calculated using 
300 Kramers pairs around ^3000. 

To summarize, we have shown that electrons in con- 
fined chaotic geometries in the presence of spin-orbit cou- 
pling exhibit spatial spin and spin current correlations 
given by the universal correlations of a Gaussian random 
spinor field. Our results can be applied to a large class of 
correlators both in the bulk and near a boundary. They 
hold not only for the correlations between amplitudes but 
also for spatial correlations of spin densities and spin cur- 
rents, more generally for any pair of local observables. A 
rigorous identity relating linear response coefficients to 
these universal correlators gives further insight into our 




FIG. 4: Spatial correlation Cj yCTz (r, r'; E) between spin cur- 
rent and spin density near a wall at x = produced, for 
instance, by applying a spin current with y polarization at 
r' — (0, 005L) and measuring the dependence with r of 
the spin polarization in z direction. Left side: numerical 
calculation from Eq. Q. Right side: universal prediction 
Cj™ M (r, r'; E) from Eq. {TO}. The parameters are the same 
as in Fig. [3] 



results and opens a straightforward possibility to mea- 
sure such correlations. 
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